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Objective. To characterize the circadian clock in 
murine cartilage tissue and identify tissue-specific clock 
target genes, and to investigate whether the circadian 
clock changes during aging or during cartilage degen- 
eration using an experimental mouse model of osteo- 
arthritis (OA). 

Methods. Cartilage explants were obtained from 
aged and young adult mice after transduction with 
the circadian clock fusion protein reporter PER2::luc, 
and real-time bioluminescence recordings were used to 
characterize the properties of the clock. Time-series 
microarrays were performed on mouse cartilage tissue 
to identify genes expressed in a circadian manner. 
Rhythmic genes were confirmed by quantitative reverse 
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transcription-polymerase chain reaction using mouse 
tissue, primary chondrocytes, and a human chondrocyte 
cell line. Experimental OA was induced in mice by 
destabilization of the medial meniscus (DMM), and 
articular cartilage samples were microdissected and 
subjected to microarray analysis. 

Results. Mouse cartilage tissue and a human 
chondrocyte cell line were found to contain intrinsic 
molecular circadian clocks. The cartilage clock could be 
reset by temperature signals, while the circadian period 
was temperature compensated. PER2::luc biolumines- 
cence demonstrated that circadian oscillations were 
significantly lower in amplitude in cartilage from aged 
mice. Time-series microarray analyses of the mouse 
tissue identified the first circadian transcriptome in 
cartilage, revealing that 615 genes (—3.9% of the ex- 
pressed genes) displayed a circadian pattern of expres- 
sion. This included genes involved in cartilage homeo- 
stasis and survival, as well as genes with potential 
importance in the pathogenesis of OA. Several clock 
genes were disrupted in the early stages of cartilage 
degeneration in the DMM mouse model of OA. 

Conclusion. These results reveal an autonomous 
circadian clock in chondrocytes that can be implicated 
in key aspects of cartilage biology and pathology. Con- 
sequently, circadian disruption (e.g., during aging) may 
compromise tissue homeostasis and increase suscepti- 
bility to joint damage or disease. 

The circadian clock governs ~24-hour cycles in 
physiology through rhythmic control of tissue-specific 
sets of clock-controlled genes (CCGs), which allows 
precise orchestration of organ function (1). In mammals 
(including humans), the circadian system is organized in 
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a hierarchical manner. The suprachiasmatic nuclei of the 
hypothalamus receive information from external time 
cues (predominantly, the light/dark cycle) and synchro- 
nize peripheral clocks in most major body organs 
through neuronal or systemic factors (2,3). The molec- 
ular basis of the circadian clock relies on the rhythmic 
activity of evolutionarily conserved clock genes and 
proteins, including those for transcriptional activators 
(Clock/Npas2 and Bmall), transcriptional repressors 
(Perl/2 and Cryl/2), and nuclear hormone receptors 
(Nrldl/2 [Rev-Erba/j3] and Rora/Rorg) (1). This 
transcriptional-translational feedback loop controls 
expression of downstream CCGs to regulate tissue- 
specific physiology. Disruption of the circadian rhythm 
(attributable to, for example, effects of aging or shift 
work) has been correlated with an increased risk of 
various human diseases, including obesity, diabetes, 
cardiovascular disease, and cancer (4). 

Articular cartilage is a specialized load-bearing 
tissue comprising an abundant extracellular matrix, 
which is produced and maintained by a sparse popula- 
tion of chondrocytes. As the only cell type residing in 
this tissue, chondrocytes have a poor capacity for endog- 
enous repair, and there is little evidence of cell division 
throughout adult life (5). Deterioration in chondrocyte 
function and survival accompanies the damage and loss 
of cartilage, a common outcome in degenerative and 
inflammatory joint diseases such as osteoarthritis (OA) 
and rheumatoid arthritis (RA). 

Several physiologic processes in cartilage exhibit 
diurnal variation, including the processes of matrix syn- 
thesis, the growth rate in the growth plate, and miner- 
alization (6-8). Mice with mutations in clock genes show 
altered regulation of bone volume (9), retarded long 
bone growth (10), and increased susceptibility to inflam- 
matory arthritis (11), providing in vivo evidence for the 
importance of the molecular clock in the skeletal sys- 
tem. Moreover, in patients with OA and those with RA, 
there is a clear circadian rhythm in the severity of pain 
and stiffness and the extent of manual dexterity (12,13). 
Expression of clock genes has been shown in bone and 
joint tissues (11,14,15) and in isolated chondrocytes 
(10,16), but no studies have directly examined the de- 
tailed mechanisms, inputs, or targets of the cartilage 
circadian clock, and none have examined how the clock 
changes during aging and in the pathogenesis of joint 
disease. 

In the present study, we demonstrate the pres- 
ence of self-sustained, functional circadian clocks in 
mouse cartilage tissue, primary chondrocytes, and a 
human chondrosarcoma cell line. Circadian transcrip- 



tome profiling revealed rhythmic expression of —3.9% 
of the genes expressed in mouse cartilage, including key 
genes associated with tissue homeostasis. In addition, 
circadian rhythms in the cartilage from aged mice were 
significantly dampened. Finally, the expression of key 
circadian clock genes was altered in articular cartilage 
from mice with experimentally induced OA. Taken 
together, our results implicate chondrocyte circadian 
clocks as a key regulator of cartilage homeostasis, and 
suggest that changes in circadian regulation may play a 
role in the pathogenesis of joint disease. 

MATERIALS AND METHODS 

Animals. All animal studies were performed in ac- 
cordance with the 1986 UK Home Office Animal Procedures 
Act. Approval was provided by local ethics review or by the 
Murdoch Children's Research Institute Animal Ethics Com- 
mittee. 

PER2::luc mice carry the firefly luciferase gene fused 
in-frame with the 3' end of the Per2 gene, creating a fusion 
protein reporter (17). Clock\19 mice harbor a deletion in exon 
19 of the Clock gene (18), producing a dominant-negative 
mutant protein. ClockA19 mice were crossed onto a PER2::luc 
background. All mice were bred in-house at the University of 
Manchester. In addition, 6-week-old male C57BL/6J wild-type 
(WT) mice were purchased from Harlan Laboratories for use 
in the time-course microarrays. 

Cartilage cultures and bioluminescence recording. 
Cartilage cultures from PER2::luc mice were prepared by 
dissection of the cartilaginous portion of the xiphoid process, 
by dissection of the femoral head cartilage from 6-10-day-old 
mice, or by scalpel microdissection of the articular cartilage 
from the femoral condyle and tibial plateau surfaces of 4.4- 
month-old mice. Cartilage was cultured on 0.4-/j,m cell culture 
inserts (Millipore), and bioluminescence was recorded in real 
time using photomultiplier tube (PMT) devices (19) or a 
LumiCycle apparatus (Actimetrics). Baseline subtraction was 
carried out using a 24-hour moving average. The circadian 
period was calculated using LumiCycle analysis software 
(Actimetrics) or using the RAP algorithm (19). Cultures were 
also visualized using a self-contained Olympus Luminoview 
LV200 microscope and recorded using a cooled Hamamatsu 
ImageEM C9100-13 EM-CCD camera. Images were obtained 
every hour for 6 days, and results were combined in ImageJ. 

For temperature entrainment studies, xiphoid tissues 
were cut into halves and cultured under PMT recorders in 
separate incubators at 37°C. After 4 days, alternating 12-hour 
square-wave temperature cycles of 38.5°C/36°C were applied; 
the temperature protocols of the 2 incubators were in anti- 
phase. After 3 full temperature cycles, cultures were returned 
to 37°C For temperature compensation studies, xiphoid cul- 
tures were incubated under PMT devices at a constant tem- 
perature of 29°C, 32°C, 37°C, and 39°C (n = 3, n = 4, n = 7, 
and n = 3, respectively). The mean period (p) was determined 
at each temperature (T). The frequency (reaction rate [R]) was 
calculated as 1/p, and the temperature coefficient (Q 10 ») was 



2336 



GOSSAN ET AL 



worked out using the following equation: Q 10 = = (Rj/Rj) ~ 10/ 
(T 2 - Tj). 

Aged animals. Male PER2::luc mice were aged to 
20-24 months (n = 10) under a 12-hour light/dark cycle. 
Young control mice were aged to —2-4 months (n = 7). 
Animals were transferred into a state of constant darkness for 
20 days. All animals were killed in the dark at their individual 
12-hour cull time (CT12) (i.e., predicted onset of activity for 
each animal). Amplitude was calculated as the peak-trough 
difference in bioluminescence of the second peak, with sub- 
traction of baseline data. A minimum of 3 peaks were used to 
determine the circadian period. Phase was determined by 
plotting the time at which the first peak occurred in culture, 
relative to CT12. For dexamethasone treatment, a final con- 
centration of 100 nM was applied directly to the recording 
medium. Amplitude was calculated from the second peak 
posttreatment. 

Histology. Xiphoid cartilage was fixed for 24 hours in 
4% paraformaldehyde, followed by decalcification in 0.8M 
EDTA (pH 7.4). Tissue was embedded in wax blocks prior to 
sectioning and hematoxylin and eosin staining using a Shandon 
Varistain 24-4 automated staining system. 

Time-series microarrays. C57BL/6J mice were en- 
trained to a 12-hour light/dark cycle before release into 
constant darkness for 39 hours prior to collection of xiphoid 
cartilage. Tissues were collected every 4 hours, commencing at 
CT3, for 40 hours. Tissue was pooled (5 animals per pooled 
sample, 3 samples per time point) and snap frozen in liquid 
nitrogen. One sample per time point was used for microarray; 
the remainder of the samples were used for quantitative 
reverse transcription-polymerase chain reaction (qRT-PCR) 
validation. 

Frozen tissue was disrupted using Mikro-Dismembrator 
S (Satorius Stedim Biotech). Extraction of RNA was carried 
out using RNeasy Fibrous Tissue Mini/Micro kits (Qiagen). 
Affymetrix Mouse430_2 GeneChips were run according to 
the manufacturer's instructions. Technical quality control 
and outlier analysis were performed using Affymetrix dChip 
software (version 2005). Background correction, quantile 
normalization, and gene expression analysis were performed 
using robust multiarray average. Raw data were deposited in 
Array Express (accession no. E-MEXP-3780). Genes with 
an average expression >2 SD higher than the mean back- 
ground of all GeneChips were classified as expressed. Circ- 
Wave Batch version 5 (provided by Dr. Roelof Hut, University 
of Groningen) was used to fit a sine-wave with 24-hour 
periodicity to each gene expression data set. The local false 
discovery rate was calculated using the fdrtool in the R 
program (results available from the corresponding author 
upon request) (20). 

Using known clock genes as a guide, a cutoff point 
of q < 0.1 was arbitrarily assigned to determine circadian 
gene expression. JTKCycle (21) was also used to identify 
circadian transcripts. A Bonferroni-adjusted P value of 0.05 
was arbitrarily set as the cutoff for significance. Functional 
annotation analysis was carried out on overlapping genes, 
using DAVID algorithms. To model the data for clustering, 
we used a Gaussian process model (22). To perform clustering, 
we specified a statistical model based on a Dirichlet process, 



a mixture of Gaussian processes (23), solved using approxi- 
mate Bayesian inference (based on the approach described in 
ref. 24). This allowed joint inference of clustering (i.e., which 
genes belong to which periodic function), alongside inference 
of the functions themselves. Validation of time-series arrays 
was carried out using TaqMan-based qRT-PCR (a list of the 
primers/probes used is shown in Supplementary Table 2, 
available on the Arthritis & Rheumatism web site at http:// 
onlinelibrary.wiley.com/doi/10.1002/art.38035/abstract), with 
results normalized to the values for p-actin expression, using 
the 2~ AAC ' method. 

Primary chondrocyte culture. Cells were prepared 
from 5-day-old PER2::luc mice as described previously (25). 
Briefly, cartilage was dissected from the femoral heads, femo- 
ral condyles, and tibial plateaus of one litter of mice and 
sequentially digested in collagenase D solution (Roche) at 
3 mg/ml and then 0.5 mg/ml. Cells were seeded directly into 
35-mm recording dishes and grown until 90-100% confluent, 
prior to synchronization and PMT recording or RNA extrac- 
tion (as described below). The chondrocyte phenotype was 
confirmed using immunofluorescence against type II collagen 
(Abeam catalog no. ab54236, with an Alexa 488-conjugated 
donkey anti-mouse antibody) and by staining of glycosamino- 
glycans in the extracellular matrix using 1% weight/volume 
Alcian blue 8GX. 

SW-1353 cell culture and lentiviral transduction. The 
SW-1353 cells used herein constitutively overexpress Sox9 to 
drive Col2al expression and maintain a chondrocyte-like phe- 
notype. For time-course experiments, confluent cells were 
synchronized with 100 nM dexamethasone for 1 hour. Samples 
were collected 24 hours postsynchronization. RNA samples 
were prepared using the RNeasy Mini kit (Qiagen). Gene 
expression was analyzed by TaqMan qRT-PCR (a list of the 
primers/probes used is shown in Supplementary Table 2, 
available on the Arthritis & Rheumatism web site at http:// 
onlinelibrary.wiley.com/doi/10.1002/art.38035/abstract), with 
results normalized to the values for GAPDH using the 2~ AAC ' 
method. For lentiviral transduction, fragments of the mouse 
Per2 (—418 bp) and Bmall (—649 bp) promoters driving 
luciferase reporters in lentiviral transfer vectors (19) were 
cotransfected with lentiviral packaging plasmids into HEK 
293FT packaging cells. Virus particles were concentrated 
in Vivaspin 20 columns (Satorius Stedim Biotech) and ap- 
plied to cells prior to synchronization and subsequent PMT 
recording. 

Induction of experimental arthritis. OA was induced 
in 10-week-old male C57BL/6 WT mice by destabilization of 
the medial meniscus (DMM) of the right knee (26), with sham 
operation of the left knee as a control. Animals were killed at 
1, 2, and 6 weeks after surgery (n = 4 animals with sham/ 
DMM-operated knees per time point). Tibial epiphyses were 
isolated, placed in RNALater (Ambion) containing 20% 
ETDA, and decalcified at 4°C for 72 hours. Total RNA was 
extracted from pooled laser-microdissected sections using 
TRIzol reagent (Invitrogen) and linearly amplified in 2 rounds 
using the MessageAmp aRNA kit (Ambion). 

Labeled RNA was hybridized to duplicate whole- 
genome microarrays (G4122A; Agilent). Raw data were pro- 
cessed in R statistical language, using the limma package. Data 
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were adjusted for multiple testing using Benjamini and Hoch- 
berg's method to control the false discovery rate. Raw data 
have been deposited in the Gene Expression Omnibus data- 
base (26). Microarrays were validated by SYBR Green qRT- 
PCR (a list of the primers used is shown in Supplementary 
Table 3, available on the Arthritis & Rheumatism web site at 
http://onlinelibrary.wiley.com/doi/10.1002/art.38035/abstract), 
using the comparative C t method. Data were normalized to the 
geometric mean expression values for 2 housekeeping genes, 
Atp5b and RpllO. 

RESULTS 

An intrinsic molecular circadian clock in car- 
tilage. To investigate circadian rhythms in cartilage tis- 
sue, we monitored the dynamics of PER2::luc protein 
expression ex vivo in xiphoid (adult mouse) cartilage 
and femoral head (6-day-old mouse) cartilage from 
PER2::luc reporter mice, using bioluminescence photon 
counting. This revealed robust circadian oscillations in 
both types of cartilage (period —24.5 hours), which 
persisted for 1-3 weeks (Figures 1A and B). (Further de- 
tails are available from the corresponding author upon 
request.) Treatment of cartilage samples (n = 3) with 
100 nM dexamethasone (a known synchronizer of peri- 
pheral clocks) immediately reinitiated a robust rhythm 
(P = 0.018 versus untreated control) (Figure 1A). 

Bioluminescence of the xiphoid tissue was also 
visualized using a high-sensitivity EM-CCD camera. This 
revealed rhythmic PER2::luc activity from the cartilagi- 
nous tissue (Figure IB) (see also Supplementary Video, 
available on the Arthritis & Rheumatism web site at 
http://onlinelibrary.wiley.com/doi/10.1002/art.38035/ 
abstract). There were no differences in the clock prop- 
erties between the femoral head cartilage and the xi- 
phoid cartilage (mean ± SD period 24.18 ± 0.3 hours 
and 23.4 ± 0.5 hours, respectively; P = 0.2 by r-test). 
Microdissected articular cartilage from the knee joints 
of adult PER2::luc mice also demonstrated circadian 
oscillations (results available from the corresponding 
author upon request). However, the small tissue size 
meant that cultures did not produce robust enough 
rhythms for further study. Therefore, xiphoid tissue was 
used as a convenient and more amenable model in 
subsequent studies. 

Because cartilage is avascular, with no direct in- 
nervation, we hypothesized that body temperature oscil- 
lations controlled by the central brain clock could be 
important in entraining local rhythms in this tissue. 
Indeed, when PER2::luc mouse cartilage was incubated 



Xiphoid 



Femoral Head 



— DexlOOnM „ -j 

HWIT I 'r~~ 

a n o -4S' 1 — . — i — . — 



4 8 

Day 



12 



0 1 2 3 4 5 6 

Day 





- 


44 hr _ 






Min ™™«™«ixj_- 


— Max 




8 10 




ii 

CT i 

if 



Before After 
lOayl) |Djy7) 



= 0.9 1 



24 £ 
1 
J8 I 



26 31 3< 41 

Temperature I C) 



Time (day) 

Figure 1. Evidence of an autonomous circadian clock in mouse 
cartilage tissue. A, Representative bioluminescence traces (in counts 
per second) in cultures of xiphoid cartilage (mean ± SD period 
24.67 ± 0.13 hours; n = 7) and articular cartilage from the femoral 
head (mean ± SD period 23.7 ± 0.6 hours; n = 4) of PER2::luc mice. 
Arrow indicates the time of treatment with 100 nM dexamethasone. B, 
Hematoxylin and eosin (H&E) staining (left) and high-sensitivity 
EM-CCD camera images (right) of xiphoid tissue. Left, A represen- 
tative H&E-stained tissue section demonstrates typical morphologic 
features of the xiphoid cartilage (first panel; bar = 500 /im). The 
second panel is a higher-magnification view of the boxed area in the 
first panel (bar = 100 jj.m). Right, EM-CCD images of the xiphoid 
tissue were obtained at 32, 44, and 56 hours post-dexamethasone 
synchronization (bars = 300 (im). C, Temperature entrainment. Two 
xiphoid cartilage cultures (represented by red and blue lines) from the 
same animal were held under antiphase temperature cycles (alternat- 
ing 12-hour cycles of 38.5°C/36°C; baseline temperature of 37°C). Left, 
Detrended trace (in counts per minute); shaded area represents 
enforcement of the temperature protocol. Right, Time of peak circa- 
dian expression on day 1 and day 7; CTO = circadian time 0 (start of 
recording). On day 7, the 2 cultures had significantly different peak 
times (P = 0.0006 by paired /-test). D, Temperature compensation. 
Left, Representative detrended traces (in counts per second) in 
cartilage (n = 4) over a 10°C temperature range. Right, Frequency 
(reaction rate), calculated as 1/period. Temperature coefficient 
(Q 10 ) = 0.9, line of best fit y = -0.00039* + 0.05422, R 2 = 0.68. 
Results are the mean ± SEM. 
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Figure 2. A less robust circadian clock in mouse cartilage during aging. A, Representative traces of PER2::luc bioluminescence from xiphoid 
cartilage of young adult mice (ages 3-4 months; n = 7) and aged mice (ages 20-24 months; n = 9). Arrow indicates time of treatment with 100 nM 
dexamethasone (Dex) (n = 4 young mice, n = 3 aged mice). B, Oscillation amplitude (left) and circadian period (right) in cartilage tissue from aged 
mice compared to young mice. * = P = 0.015 by f-test. NS = not significant. C, Circadian phase at the time between 12-hour cull and the first peak 
of bioluminescence in culture. Vertical lines represent the mean for each group (P = 0.46). D, Similar extent of induction of oscillation amplitude 
by dexamethasone treatment in young and aged mouse cartilage tissue; amplitude is expressed as a percentage of the amplitude of the young tissue 
before treatment. Results are the mean ± SEM. Data were analyzed by two-way analysis of variance. There were no significant effects of age group 
(P = 0.77), and interactions were not significantly different (P = 0.9992). ** = P = 0.003. 



under alternating 12-hour temperature cycles of 38.5°C/ 
36°C, oscillation amplitude was induced (Figure 1C). 
Incubating cultures of cartilage tissue from the same 
animal under antiphase temperature cycles drove the 
tissue into antiphase (Figure 1C). Altered phase was 
maintained under subsequent conditions of constant 
temperature, thus demonstrating that the observed ef- 
fect on the circadian rhythm was a clock-based change 
rather than a passive response to changes in tempera- 
ture. 

Recent evidence suggests that temperature fluc- 
tuations will target circadian clock genes at the tran- 
scriptional level (3,27). To this end, we lentivirally 
transduced a human chondrosarcoma cell line with clock 
promoter luciferase reporters, and then administered 
single 6-hour pulses of heat (+2.5°C) or cool (— 2.5°C). 



Consistently, we found that the transcriptional levels of 
Per2::\\ic and Bmallr.hic were altered (results available 
from the corresponding author upon request). 

A defining feature of a circadian oscillator is 
temperature compensation, the maintenance of robust 
periodicity over a broad temperature range (28). Indeed, 
our data demonstrated that the circadian rhythm in 
cartilage was effectively temperature compensated. The 
circadian period ranged from 23.5 hours to 26.5 hours 
over a 10°C temperature range (temperature coefficient 
Q 10 = 0.9) (Figure ID). Taken together, these findings 
establish the presence of a temperature-compensated, 
autonomous circadian clock in cartilage tissue and sug- 
gest that this clock senses suprachiasmatic nuclei-driven 
daily changes in body temperature as one means of 
maintaining synchrony with the rest of the body. 
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Figure 3. Identification of the circadian transcriptome in mouse cartilage. A, Heat map depicting the expression level of the 615 circadian genes 
(represented by 705 probe sets), showing that 3.9% of the cartilage transcriptome is circadian, as identified using CircWave Batch (local false 
discovery rate < 0.1, equivalent to approximate R 2 > 0.7) and JTKCycle (Bonferroni-adjusted P < 0.05). Genes are organized according to timing 
of peak expression. B, Clustering and identification of periodic functions. Each panel represents a single cluster, with normalized log-expression 
levels marked as connected points for each gene over time. Colored lines represent the inferred cluster-wise periodic function, with 2 standard 
deviations marked by the shaded area. Clusters fall into 4 broad categories, based on the time of peak expression. In A and B, gray-shaded bars 
represent subjective day, and solid black bars represent subjective night. C, Tissue specificity of cartilage clock-controlled genes. Venn diagram shows 
a comparison of cartilage circadian genes with previously published gene lists from studies of the liver and skeletal muscle of mice. The total number 
of genes identified as circadian in each tissue is shown in parentheses; areas of overlap indicate common genes. 



Impact of aging on the cartilage circadian 
rhythm. With aging, systemic circadian rhythms in body 
temperature and hormone release are altered (29). This 
may consequently alter the cartilage clock. To test this, 
we compared PER2::luc rhythms in cartilage from aged 
mice (ages 20-24 months) and young adult mice (ages 
2-4 months). Although there were no significant changes 
in the circadian period or phase, oscillation amplitude 
was significantly reduced (reduction of —40%) in aged 
mouse cartilage (Figures 2A-C). Interestingly, upon 
dexamethasone treatment, robust rhythmicity was re- 
initiated to a similar extent (increase of ~5-10-fold) 
in both aged and young adult mouse cartilage tissue 



(Figure 2D), indicating that intrinsic pacemaking mech- 
anisms within aged cartilage may still be intact. These 
studies provide evidence that the circadian rhythm in 
cartilage loses its robustness with advanced age. 

Identification of the circadian transcriptome in 
cartilage. To gain insight into the function of the circa- 
dian clock in cartilage, we aimed to identify its tissue- 
specific targets. To this end, we performed a time-series 
microarray using xiphoid tissue from mice kept in con- 
stant darkness. This eliminated the effect of the light/ 
dark cycle, and allowed us to reveal endogenously 
circadian genes. Time points were every 4 hours for 
40 hours. To identify circadian transcripts, we used 2 
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well-recognized algorithms based on different statistical 
models. With use of the CircWave Batch algorithm, we 
identified 896 genes oscillating with a period of —24 
hours (q < 0.1). Parallel analysis was performed using 
the JTKCycle method (21), which identified 1,574 rhyth- 
mic genes with a 20-28-hour period (adjusted P < 0.05). 
To be stringent, we considered only those genes identi- 
fied by both methods as circadian. Thus, we identified 
615 genes with rhythmic expression and a circadian 
period of —24 hours (3.9% of expressed genes, repre- 
sented by 704 probes) (Figure 3A and Supplementary 
Table 1, available on the Arthritis & Rheumatism web site 
at http://onlinelibrary.wiley.com/doi/10.1002/art.38035/ 
abstract). 

Rhythmic cartilage genes were clustered accord- 
ing to their inferred periodic functions, determined 
using a Gaussian process model. Genes fell into 11 
distinct clusters, which were grouped into 4 categories 
based on the phase of their peak expression (Figure 3B). 
Most genes peaked during subjective daytime (Figure 
3A and results available from the corresponding author 
upon request). Functional annotation and gene ontology 
(GO)-term analysis using DAVID revealed several 
overrepresented functional groups (a full list is available 
from the corresponding author upon request), with the 
top overrepresented groups being "purine/ATP bind- 
ing," "circadian rhythm," "extracellular matrix location/ 
microfibril," "negative regulator of (macromolecule) 
biosynthetic processes," "proteolysis," and "apoptosis." 
Interestingly, 81.3% of genes related to "negative regu- 
lation of biosynthesis" and 78.6% of genes related to 
"proteolysis" peaked during the animals' resting phase 
(results available from the corresponding author upon 
request). 

We compared our gene list to published sequen- 
tial chromatin immunoprecipitation assay data on the 
expression of circadian clock transcription factor Bmall 
in the liver (30), and we found that 132 of our rhythmic 
genes had previously determined Bmall binding sites 
(see Supplementary Table 1, available on the, Arthritis & 
Rheumatism web site at http://onlinelibrary.wiley.com/ 
doi/10.1002/art.38035/abstract). Only 19 circadian genes 
were common (Figure 3C) between our gene expression 
data set from mouse cartilage and previously published 
data from the liver (31) and skeletal muscle (32) of mice, 
most of which were core clock genes, supporting the 
concept of tissue-specific clock function. 

The expression profiles of canonical clock genes 
(Bmall, Per2, Cryl, Nrldl, Dbp, and E4bp4) and se- 
lected cartilage genes (Adamts4, Adamts9, and Mmpl4) 
were validated by temporal qRT-PCR in mouse xiphoid 
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Figure 4. Validation of circadian genes in mouse cartilage. A, Rela- 
tive expression of a panel of clock genes and putative clock-controlled 
genes in mouse xiphoid tissue, as determined by microarray and 
validated by quantitative reverse transcription-polymerase chain reac- 
tion (qRT-PCR). Data were normalized to the nadir of expression. All 
genes tested followed the expression pattern identified by the array; 
when analyzed using one-way analysis of variance, all genes showed 
significant change with time (P < 0.0001; for MrnpU, P = 0.0039). 
Gray-shaded bars denote subjective day, and solid black bars denote 
subjective night. Results are the mean ± SEM (n = 3 replicates, 
pooled from 5 animals per replicate). B, Abolition of circadian 
rhythmicity in Clock A19 mutant animals. Left, Xiphoid tissue cultures 
from ClockA19 mice on a PER2::luc background showed no molecular 
rhythmicity, as compared to wild-type (WT) mice, in cultures (results 
representative of 2 independent experiments; bioluminescence ex- 
pressed in counts per second). Right, Time-dependent expression of 
clock output genes was abolished in Clockkl9 animals. The qRT-PCR 
data were normalized to the nadir of expression in WT mice. Results 
are the mean ± SEM (n = 5). ** = P < 0.01; ■**■*= P < 0.001, by 
Mest. NS = not significant; ZT = zeitgeber time; ZTO = lights-on. 



tissue (Figure 4A). To confirm the importance of the 
molecular clock in controlling cartilage gene rhythmic- 
ity, we used ClockA19 mutant mice on a PER2::luc 
background. The expression and behavior of genes in 
ClockA19 mice are known to be arrhythmic (18). We 
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Figure 5. Cell-autonomous circadian clocks in chondrocytes. A, Representative circadian oscillations of bioluminescence from immature articular 
chondrocytes isolated from 5-day-old PER2::luc mice (n = 5). B, Confirmation of the chondrocyte phenotype. Left, Cell morphology under phase 
contrast. Bar = 100 /j.m. Middle, Alcian blue staining of glycosaminoglycans in the extracellular matrix. Bar = 50 jum. Right, Immunofluorescence 
for type II collagen (green) with DAPI nuclear counterstain (blue). Bar = 50 jum. C, Time-dependent expression of clock genes in primary 
chondrocytes. Cells were synchronized and gene expression was assayed by quantitative reverse transcription-polymerase chain reaction (qRT-PCR) 
at 6-hour intervals. Data were normalized to the nadir of expression. Results are the mean ± SEM (n = 3 technical replicates). P < 0.001 for all 
genes, by one-way analysis of variance (ANOVA). D, Expression of Per2::\uc and Bmallvluc reporters in separate cultures of SW-1353 cells after 
lentiviral delivery. Cells were synchronized by treatment with dexamethasone for 1 hour before bioluminescence recording. Representative 
detrended traces are shown (mean ± SEM period 24.1 ± 0.3 hours; n = 4). E, Time-dependent expression of cartilage clock targets in human 
chondrosarcoma cells. Cells were synchronized, and gene expression was assayed by qRT-PCR at 4-hour intervals. Data were normalized to the nadir 
of expression. Results are the mean ± SEM (n = 4). P < 0.05 for all genes, by one-way ANOVA. 



established that PER2::luc rhythmicity was lost in the 
mutant mouse cartilage, and that time-dependent ex- 
pression of the circadian output genes Dbp and E4bp4 
was abolished (Figure 4B). Importantly, time-dependent 
expression of the cartilage genes Adamts4 and Mmpl4 
was similarly abolished in the ClockA19 mouse cartilage 
(Figure 4B), supporting the conclusion that these genes 
are clock-controlled. 

Cell-autonomous circadian rhythms in chon- 
drocytes. To confirm that chondrocytes harbor cell- 
autonomous circadian clocks, we isolated immature 
articular chondrocytes from 5-day-old PER2::luc mice. 
PMT recording revealed self-sustained circadian oscilla- 



tions resembling those from the cartilage tissue (Figure 
5A). The mouse cells demonstrated typical morpho- 
logic features of chondrocytes, characterized by Alcian 
blue staining of glycosaminoglycans in the extracellular 
matrix and high Co/2 expression (Figure 5B). Further- 
more, clock-synchronized primary chondrocytes time- 
dependently expressed the core clock genes Nrldl, Per2, 
Bmall, and Cryl (Figure 5C). 

To test whether human chondrocytes also exhibit 
circadian rhythmicity, we used lentiviral delivery to 
stably transduce a human chondrosarcoma cell line 
(SW-1353) with lucif erase reporters driven by promoter 
fragments of the clock genes Per2 or Bmall. Both 
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reporters oscillated robustly over several days (Figure 
5D), in antiphase to one another. Similar to that ob- 
served in mouse cells, the clock genes Nrldl, Per2, and 
Bmall were time-dependently expressed in SW-1353 
cells (Figure 5E). Interestingly, putative clock target 
genes (Adamts4/9, Mmpl4, and Gnl3) were also ex- 
pressed in a time-dependent manner (Figure 5E), al- 
though Adamts4 expression was less robust. 

Taken together, these data support the concept 
that chondrocytes contain a cell-autonomous molecular 
clock capable of driving downstream targets, extending 
our findings from a mouse model into human cells. This 
implies that there is a possible conservation of clock 
mechanisms in chondrocytes from different species. 

Clock gene changes in articular cartilage in an 
experimental mouse model of OA. We next sought to 
address whether joint articular cartilage expresses clock 
genes, and whether the expression of clock genes is 
altered in a surgical DMM model of OA. Laser micro- 
dissection allowed us to focus on lesion sites, to study 
factors initiating damage. One week post-DMM, carti- 
lage tissue damage was minimal; 2 weeks post-DMM, 
tissue exhibited early focal degeneration; 6 weeks 
post-DMM, tissue was more extensively damaged, with 
fibrillation. 

Our microarray data set of the gene networks 
involved in mouse OA is presented elsewhere (26). We 
interrogated these data and found that the canonical 
clock genes Bmall, E4bp4, Nrldl, Per2, Cryl, Dbp, and 
Rora were indeed expressed in the DMM articular 
cartilage. Bmall was expressed at a significantly higher 
level in DMM cartilage than in sham control cartilage 
at 1, 2, and 6 weeks post-DMM (P < 0.05) (Figure 6). 
The clock output genes Nrldl, Dbp, Rora, and E4bp4 
were also significantly altered 2 weeks post-DMM (P < 
0.05) (Figure 6), consistent with the time point at which 
the greatest change in Bmall expression was observed. 
Changes in clock gene expression were validated using 
qRT-PCR (results available from the corresponding 
author upon request). Taken together, these data indi- 
cate that clock gene expression changes early in the 
progression of experimental OA. 

DISCUSSION 

Mice with mutations in clock genes show ec- 
topic ossification of tendon and cartilage, altered regu- 
lation of bone volume, an increased susceptibility to 
inflammatory arthritis, and, as recently reported, de- 
fects in long bone growth through clock control of 
chondrocyte differentiation (9-11,33). These studies 
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Figure 6. Deregulation of clock genes in an experimental mouse 
model of osteoarthritis. Articular cartilage mRNA expression of core 
clock and clock output genes was determined by microarray analysis. 
Results of all probes for each gene are shown, with findings expressed 
as the fold up-regulation or down-regulation of gene expression in 
joints subjected to destabilization of the medial meniscus (DMM) 
compared to sham-operated joints in mice at 1, 2, or 6 weeks 
post-DMM surgery. Values are the mean fold change in 4 mice. The 
shaded region between + 1-fold and —1-fold represents no differential 
expression. * = adjusted P < 0.05; ** = adjusted P < 0.01. 

highlight the potential importance of circadian clocks in 
the skeletal system. However, to the best of our knowl- 
edge, no studies have demonstrated functional and 
cell-autonomous circadian clocks in fully differentiated 
cartilage tissue. Therefore, we set out to address the 
following questions. 1) Do cartilage and chondrocytes 
contain self-sustained circadian clocks? 2) What factors 
entrain the cartilage clock? 3) What cartilage-specific 
genes/pathways are controlled by the circadian clock? 
4) How does the cartilage circadian clock change 
with aging and during the pathogenesis of experimental 
arthritis? 

Herein, we have demonstrated cell-autonomous 
circadian clocks in cartilage tissue and chondrocytes. 
Similar to previously described peripheral circadian 
clocks (3,28), the cartilage clock is entrained by temper- 
ature cycles approximating body temperature fluctua- 
tions and its periodicity is effectively maintained across 
a wide temperature range. Thus, sensing changes in body 
temperature may be one way for the cartilage clock to 
maintain synchronization to the rest of the body. Al- 
though chondrocytes are dispersed in the avascular and 
aneural cartilage, systemic factors freely diffuse and it is 
probable that multiple entraining factors, such as hor- 
mone signaling, act alongside temperature to entrain the 
clock. 

Age is the major risk factor for developing OA. 
Its precise role in disease pathogenesis is unknown. 
Disruption to whole-body circadian rhythms through 
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aging could alter systemic time cues and compromise 
cartilage clock function. Indeed, the results of the pres- 
ent study show that PER2::luc circadian rhythms were 
dampened in aged mouse cartilage tissue. We speculate 
that age-associated dampening in circadian rhythms 
could disrupt downstream OA-associated genes and 
increase disease susceptibility. Circadian amplitude in 
aged tissue could be restored to that of young cartilage 
by treatment with dexamethasone, suggesting that in- 
herent changes in the molecular clock are unlikely. 
Instead, we hypothesize that age-related dampening of 
the cartilage clock may be a consequence of altered 
systemic time cues. This is consistent with a previous 
report showing no difference in the molecular clock of 
cultured human fibroblasts between old and young 
individuals, although their serum had different clock- 
synchronizing properties (34). Our data suggest that 
the aged cartilage clock might be re-tuned by factors 
strengthening systemic circadian rhythms, such as sched- 
uled exercise (35), restricted meal times (36), or the 
scheduled warming and cooling of joints; this may have 
significant impact on the future management of joint 
diseases such as OA. 

Peripheral circadian clocks temporally regulate 
tissue-specific functions through control of downstream 
CCGs. Indeed, 3-15% of total transcripts are expressed 
rhythmically, depending on the tissue, sampling fre- 
quency, and analysis stringency (31,37). Herein, we have 
identified the first circadian transcriptome profile in 
cartilage tissue. Our results revealed rhythmic expres- 
sion of —3.9% of all expressed genes, with the majority 
of the genes being tissue specific. By selecting genes 
identified by both CircWave and JTKCycle, we gener- 
ated a relatively conservative circadian gene list. How- 
ever, it is likely that some of the genes identified are 
directly responsive to time cues, such as body tempera- 
ture, rather than being local clock driven. Studies in 
other tissues (38,39) have demonstrated only a small 
number of such system-driven genes (10-16%). Indeed, 
we demonstrated cell-autonomous oscillations of clock 
genes and putative targets in cells lacking systemic time 
cues. 

The rhythmic transcripts identified in cartilage 
exhibited 11 distinct patterns of expression. We predict 
that genes within each cluster are coregulated by similar 
clock proteins and tissue-specific transcription factors. 
Interestingly, in our mouse cartilage data set, early-day 
transcripts included several extracellular matrix com- 
ponents (fibrillins, laminin, collagens, and netrin), and 
almost one-half of the rhythmic proteolysis genes 
peaked at this early-day time point. Restricting genes 



controlling macromolecule turnover to the resting phase 
may confer a selective advantage, allowing more effec- 
tive tissue remodeling and repair following (nocturnal) 
activity bouts. It is possible that disruption to this 
time-coordinated homeostasis during aging increases 
OA susceptibility. 

Disruptions in several of the rhythmic genes/ 
pathways identified herein have previously been linked 
to joint diseases (40-44). These include genes involved 
in catabolic signaling (Adamts4, Adamts9, and Mmpl4), 
anabolic signaling (Tgff3r3 and Timp4), and the oxidative 
defense pathway (Gclc and Gstml), as well as apoptosis- 
related genes (Xiap). Some genes have even been sug- 
gested as drug targets, for example, the Adamts family 
of aggrecanases (45). Interestingly, a newly identified 
OA susceptibility gene, Gnl3 (42), was also rhythmic. 
The circadian nature of these genes may help inform 
appropriate targeting strategies, including, for example, 
timed drug administration (chronotherapy), with the 
aim of improving efficacy and reducing toxicity (46). 
Our study also raises the possibility of targeting the gene 
activity of chondrocytes through the use of circadian 
clock-acting compounds (47,48), which may prove of 
use in the treatment of diseases such as OA. 

It should be noted that we used xiphoid tissue as 
a model for mature cartilage in the time series, due to 
limitations involved in microdissecting articular carti- 
lage, but that both tissues exhibit circadian rhythms 
(results available from the corresponding author upon 
request). Furthermore, we demonstrated that both mouse 
primary chondrocytes and human chondrosarcoma- 
derived SW-1353 cells rhythmically express clock genes 
in vitro (Figure 5), although we recognize that the lat- 
ter cell line is a limited model for primary cells or 
tissue (49). 

Results of previous studies have suggested that 
clock genes change during inflammatory joint disease 
(11,15), but no studies have addressed whether clock 
genes in articular cartilage change during OA initiation/ 
progression. Herein, we have shown that clock genes are 
indeed expressed in articular cartilage. Moreover, dur- 
ing the development of OA in a mouse DMM model, the 
expression of these genes was disrupted. Interestingly, 
most of the changes were seen during the early stage of 
the disease (2 weeks post-DMM), not when the disease 
phenotype was more established (6 weeks post-DMM). 
This could mean that clock gene changes contribute to 
the disease phenotype through the disruption of CCGs, 
rather than being a consequence of tissue damage. Due 
to technical limitations, we compared expression 
changes of clock genes at a single daytime time point 
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(10:00 AM-noon). Future studies with more time points 
will be required to determine changes in the period or 
amplitude of these genes. 

Taken together, our data provide substantial 
insight into the circadian regulation of cartilage biologic 
processes. These results provide a firm base upon which 
to further explore the link between the circadian clock, 
aging, cartilage homeostasis, and OA. 
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